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Abstract 

This paper reports the finding of a simple one-parameter family of three-dimensional quadratic au- 
* tonomous chaotic systems. By tuning the only parameter, this system can continuously generate a variety 

of cascading Lorenz-like attractors, which appears to be richer than the unified chaotic system that con- 
' ' tains the Lorenz and the Chen systems as its two extremes. Although this new family of chaotic systems 

has very rich and complex dynamics, it has a very simple algebraic structure with only two quadratic 
terms (same as the Lorenz and the Chen systems) and all nonzero coefficients in the linear part being 
— 1 except one —0.1 (thus, simpler than the Lorenz and Chen systems). Surprisingly, although this new 
system belongs to the family of Lorenz-type systems in some existing classifications such as the general- 
ized Lorenz canonical form, it can generate not only Lorenz-like attractors but also Chen-like attractors. 
This suggests that there may exist some other unknown yet more essential algebraic characteristics for 
* classifying general three-dimensional quadratic autonomous chaotic systems. 

o 

^ 1 Introduction 

• • 

The scientific story of chaos dates back to earlier 1960s [l], when Lorenz studied the atmospheric convection 
phenomenon he found a chaotic attractor in a three-dimensional (3D) quadratic autonomous system [2 . The 
now-classic Lorenz system is described by 

d 

X = a {y — x) 

y = rx — y — xz (1) 
z = —bz + xy, 

which is chaotic when a = 10, r = 28, ^ = §• 

Ever since then, many researchers were wondering and pondering whether the discovery of the Lorenz 
system was just a lucky incident or there actually exist other closely related systems around it? 
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In 1999, from an engineering feedback anti-control approach, Chen provided a certain answer to this 
question by his finding of a new system [3l[4], lately referred to as the Chen system by others, described by 

X = a{y — x) 

y = {c — a)x — xz -\- cy (2) 
i = —bz + xy^ 

which is chaotic when a = 35, 6 = 3, c = 28. Whereafter, it has been proved that Chen's attractor exists 
[5] and the Chen system displays even more sophisticated dynamical behaviors than the Lorenz system [6 . 

It is also interesting to recall a unified chaotic system, which was constructed to encompass both the 
Lorenz system and the Chen system [7 . This unified chaotic system is by nature a convex combination of 
the two systems, and is described by 

x = (25a + 10)(7/-x) 

y= {28- 35a)x - xz ^ (29a - l)y (3) 
z — —^^z + xy. 

When a = it is the Lorenz system while with a = 1 it is the Chen system, and moreover for any a G (0, 1) 
the system remains to be chaotic. 

As a result of several years of continued research endeavor along the same line, a family of generalized 
Lorenz systems was found and characterized O [H [TOl [11] , which is defined through the so-called generalized 
Lorenz canonical form, as follows: 

X \ / an ai2 ^ \ f ^ \ ( ^ \ 

y \ = \ a2i a22 ^ + ' (4) 

z ) VO ass J \ z J \ xy ) 

where x^ y^ z are the state variables of the system. This canonical form contains a family of chaotic systems 
with the same nonlinear terms, the same symmetry about the z-axis, the same stability of three equilibria, 
and similar attractors in shape with a two-scroll and butterfly-like structure 

According to [8 (see also [12 ), the Lorenz system satisfies ai2a2i > while the Chen system was found 
to satisfy ai2a2i < 0. In this sense, the Chen system is dual to the Lorenz system. In between, it was also 
found that there is a transition, the Lii system, which satisfies ai2a2i = [13 . 

Lately, another form of a unified Lorenz-type system and its canonical form were developed, which contain 
some generalized Lorenz-type systems and their corresponding conjugate Lorenz-type systems [HI [15]. 

Moreover, another classification was developed in [16 , where system (3) is classified into two groups, the 
Lorenz system group if aiia22 > and the Chen system group if aiia22 < 0. This classification leads to the 
finding of a transition, the Yang-Chen chaotic system, which satisfies aiia22 = 0. 

Comparing these two classification, one can see that every such system is classified according to its 
algebraic structure. Specifically, it is determined either by ai2a2i or by aiia22. 

Given the above discussions, at this point it is interesting to ask the following questions: 

1. Concerning the conditions on the signs of ai2a2i, defined in [8 (see also [12 ), are these signs essential 
to the system dynamics? For instance, can a system belonging to the Lorenz systems family generate 
Chen-like attractors? 

2. Concerning the conditions on the signs of aiia22, defined in [16 , are these signs essential to the system 
dynamics? For instance, can a system belonging to the Chen system group generate Lorenz-like 
attractors? 
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3. Concerning the unified chaotic system ([3|, is there any other simpler chaotic system that can also 
generate both Lorenz-like attractor and Chen-like attractor? 



This paper aims to provide certain answers to the above questions via the finding of a new family of 
Lorenz-like systems. This new family has only one real parameter in a simple algebraic form but demonstrates 
very rich and complex dynamics in a way similar to the generalized Lorenz canonical form [9j [lOl [11] . It 
belongs to the Lorenz-type of systems defined in [8] and it is also classified into the Chen system group 
defined in [16 . However, it has a very simple form and can generate both visually Lorenz-like attractors 
and Chen-like attractors by gradually changing the single parameter. This reveals that further study of 
the system algebraic structure and its effects on the system dynamics remains an important and interesting 
challenge. 

2 The new family of chaotic systems 

The new system under investigation is described by 



where r is a real parameter. 

This is a one-parameter family of chaotic systems in the sense that as the real parameter r is gradually 
varied, a sequence of chaotic attractors can be continuously generated from the system, with very rich and 
complicated dynamics, as further demonstrated below. 

Before proceeding to chaotic dynamics, some basic dynamical properties of this new system is analyzed 
following a standard routine, which is necessary for understanding the nature of the new system. 

2.1 Symmetry and dissipativity 

First, it is apparent that system ([5| has a natural symmetry about the z-axis under the following transfor- 
mation: 




(5) 



S{x,y,z) {-x,-y,z). 
Second, one may construct the following Lyapunov function: 




(6) 



which gives 
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the system is dissipative under the condition of r < 1.1. More precisely, since V{t) = V{to)e ^^'^ any 
initial volume Vq containing the system trajectories shrinks to zero as t ^ +00 at an exponential rate of 
1.1 - r. 



2.2 Equilibria and their stability 

It is obvious that the origin is a trivial equilibrium of system ([5|. Other nonzero equilibria can be found by 
solving the following equations simultaneously: 



-x-y = 0; 



-X -\- ry — xz = 0; 



When r > — 1, system (5) has three equilibria: 0(0, 0, 0), £^1 ^ 
By linearizing system (5) at 0(0,0,0), one obtains the Jacobian 



OAz^xy = 0. 

10 ' ^ ' y ' ^2 \^ Y 10 ' Y 10 ' ^ 



J\o = 



-1 

-I- z 

y 



r 

X 





— X 

-0.1 



(8) 



whose characteristic equation is 

det(A/ - J\o) = + (1.1 - r)A^ - (l.lr + 0.9)A - 0.1(r + 1) = 0, 

which yields 



Ai = -0.1 < 0, 
A2 = -0.5 + 0.5r 

A3 = -0.5 + 0.5r 



V5H 


h 2r + 


2 


V5H 


h 2r + 



>o, 

<0. 



Obviously, the equilibrium 0(0, 0, 0) is a saddle point, at which the stable manifold W^{0) is two-dimensional 
and the unstable manifold W^(0) is one-dimensional. 

Similarly, linearizing the system with respect to the other equilibria, Ei and £^2, yields the following 
characteristic equations: 

det(A/ - J\e,,^) = A^ + (1.1 - r)X^ + 0.2A + 0.2(r + 1) = 0. 

According to the Routh-Hurwitz criterion, the following constraints are imposed: 

Ai 



A2 = 
A3 = 



(1.1 -r) > 0, 

1.1 -r 0.2(r+l) 
1 0.2 

0.2(r + l)A2 >0. 



= 0.1 - 2r > 0, 



(9) 
(10) 

(11) 



This characteristic polynomial has three roots, all with negative real parts, under the condition of — 1 < r < 
0.05. Therefore, the two equilibria Ei and E2 are both stable nodes, or node-foci, if — 1 < r < 0.05. 

However, if the condition 0.05 < r < 1.1 holds, then the characteristic equation has one negative real 
root and one pair of complex conjugate roots with a positive real part. Therefore, the two equilibria Ei^2 
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are both saddle- foci, at each of which the stable manifold W^{Ei^2) is one-dimensional and the unstable 
manifold VK^(^i,2) is two-dimensional. 

Summarizing the above analysis and discussions establishes the following result: 

Theorem 2.1. With parameter —1 < r < 1.1, system has three equilibria: 



0(0,0,0), 




Furthermore, 

(i) 0(0,0,0) is a saddle point, at which the stable manifold W^{0) is two-dimensional and the unstable 
manifold W^(0) is one- dimensional; 

(ii) if 0.05 < r < 1.1, then the two equilibria Ei and E2 are both saddle-foci, at each of which the stable 
manifold W^(£^i,2) is one- dimensional and the unstable manifold W^(£^i,2) is two-dimensional; 

(iii) if —1 < r < 0.05, then the equilibria Ei and E2 are both stable nodes, or node-foci, at each of which the 
stable manifold W^{Ei^2) is three-dimensional. 

The Jacobian eigenvalues evaluated at the three equilibria of the unified chaotic system (|3| and of the 
new system ([5| are listed in Table 1 for comparison. 

Given the above analysis, system ^ will be discussed only within the parameter region of 0.05 < r < 1.1 
below. 

2.3 Remarks on the classification of the new system 

On one hand, it is noted that in system ([5| one has ai2a2i = 1 > 0, so it belongs to the Lorenz system 
family defined in [8]. 

On the other hand, by making the following transformation: 

T : {x,y,z) {x,-y,-z), 

system ^ becomes 

{X = —X + y 
y = X -\- ry — xz (12) 
i = —O.lz + xy^ 

which satisfies aiia22 < and an < 0, thus belong to the Chen system group according to the classification 
in [16] (see Tables 3 and 4 therein). 

It is therefore very interesting to see that the new system ([5| belongs to the same class (either Lorenz-type 
of systems defined in [8] or Chen system group defined in [16 ), yet can generate both Lorenz-like attractors 
and Chen-like attractors by gradually changing the single parameter r from 0.05 to 0.74, as further detailed 
next. 



3 Chaotic behaviors and other complex dynamics 

A complete transition from Lorenz-like to Chen-like attractors in system ([5| is shown in Figjl] The bifurca- 
tion diagram with respect to parameter r is shown in Figj2] From these figures, one can observe that system 
([5]) evolves to a chaotic state and eventually to a limit cycle, as the parameter r gradually increases. 



5 



Figure 1: Colored figures showing the projected orbit transition from Lorenz-hke to Chen-hke attractors in 
system ([5|. 



1.5 




-1 - 



-0.2 □ 0.2 0.4 0.6 0.8 



Figure 2: Bifurcation diagram of system ([5| with r G [—0.2,0.8]. 

To verify the chaotic behaviors of system ([5|, its three Lyapunov exponents Li > L2 > Ls and the 
Lyapunov dimension are calculated, where the latter is defined by 

in which j is the largest integer satisfying Xli=i ^ X]i=i < 0. 

Note that system ([5| is chaotic if Li > 0, L2 = 0, L3 < with \Li\ < \Ls\. Figjs] shows the dependence 
of the largest Lyapunov exponent on the parameter r. In particular, for several values of r, the Lyapunov 
exponents and dimensions of the system ([5| are summarized in Table 2. From Figjsj it is clear that the 
largest Lyapunov exponent increases as the parameter r increases from 0.05 to 0.7. 

Case 1: r = 0.05 

When r < 0.05, the corresponding system has one saddle and two stable node- foci (for reference, see [16 ). 
The case of r < 0.05 is subtle and somewhat complicated, which will be further studied elsewhere in the 
near future. 

Here, the discussion starts from r = 0.05. When r = 0.05, this new system has one saddle and two 
centers. The two centers have the same set of eigenvalues: one negative real eigenvalue and two conjugate 
complex eigenvalues with zero real parts. Nevertheless, the Lyapunov exponents are Li = 0.03467, L2 = 0, 
and L3 = —1.0843, and the Lyapunov dimension is Dl = 2.0317, for initial values (1, 1, 1). This convincingly 
implies that the system is chaotic. 
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0.05 0.22 0.39 0.5G 0.73 



Figure 3: The largest Lyapunov exponent of system ([5| with r G [0.05,0.75]. 
Case 2: r G (0.05,0.7] 

In this case, the new system has one saddle and two saddle-foci. The two saddle-foci have one negative real 
eigenvalue and two conjugate complex eigenvalues with a positive real part. Moreover, the largest Lyapunov 
exponent is larger than zero, which is increased as the parameter r is increased from 0.05 to 0.7. At the same 
time, the attractor of this system is changed from Lorenz-like to Chen-like as the parameter r is increased 
from 0.05 to around 0.7, as seen in Figj4]^a)-5(f). 

Case 3: r G (0.7,0.74) 

When r lies in this interval, system ([5| is still chaotic, for Figjs] shows that the largest Lyapunov exponent 
is positive. However, the attractor is neither Lorenz-like nor Chen-like, as shown in FigQg), when r = 0.72. 
Moreover, the chaotic property becomes less significant as compared with Case 2. Indeed, the largest 
Lyapunov exponent decreases when the parameter is increased from r = 0.7 to r = 0.74, as shown in Figj3| 

Case 4: r G [0.74, 1] 

In this case, although the new system has one saddle and two saddle-foci, associated with one negative 
real eigenvalue and two conjugate complex eigenvalues with a positive real part, the system is not chaotic. 
Instead, the system has a limit cycle, though different from Case 2 and Case 3. Fig|4jh) displays the limit 
cycle when r = 0.75. 

4 Comparison between the new family and the unified system 

It is interesting to compare the new family of chaotic systems with the now-familiar generalized Lorenz 
systems family, particularly the corresponding Lorenz- type and Chen- type systems, as shown in Figs j6] and 

13 

Recall in particular the unified chaotic system (|3|, which is chaotic for all a G [0,1], with a positive 
largest Lyapunov exponent as shown in Figjs] 
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Figure 4: Attractors of system ([5| with parameter values: (a) r = 0.2; (b) r = 0.3; (c) r = 0.4; (d) r = 0.50; 
(e) r = 0.6; (f) r = 0.7; (g) r = 0.72; (h) r = 0.75. 




Figure 5: The largest Lyapunov exponent of the unified chaotic system with a G [0, 1]. 



There are several common features between the new family of chaotic systems ([5| and the unified chaotic 
system ([3|: 

1. They are both one-parameter family of chaotic systems, representing infinitely many non-equivalent 
chaotic systems. 

2. They both demonstrate a continuously changing process generating from Lorenz-like to Chen-like 
attractors. 



3. They both have a similar simple algebraic structure with two same nonlinear terms (i.e., - 
second equation and xy in the third equation), juts like the Lorenz and the Chen systems. 



the 



4. They both have the same 2;-axis rotational symmetry, therefore all their attractors are visually similar 
in shape with a two-scroll and butterfiy-like structure. 

5. They both have three equilibria: one saddle and two saddle-foci. Moreover, the two saddle- foci have 
the same eigenvalues: one negative real number and two conjugate complex numbers with a positive 
real part. Specifically, for the unified chaotic system, the positive real part of the conjugate complex 
eigenvalues increases from the Lorenz system (0.0940) to the Chen system (4.2140), while for the new 
system, it increases as r is increased from 0.05 to 0.7. 

The main difference between these two one-parameter families of chaotic systems are quite subtle. For 
system ([5|, the real part of its conjugate complex eigenvalues is precisely zero when r = 0.05. This is a 
critical point when the stability of the equilibrium pairs is changed from stable to unstable. Therefore, one 
may define the value of r = 0.05 as a boundary, or more precisely as a starting point of the Lorenz systems 
family. Since the parameter r can take values on both sides of r = 0.05 to generate chaos, the new family 
system ([5| is richer and more interesting than the unified system system (|3|. 
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(g) 11 (h) 



Figure 6: Attractor of system ([5| with r = 0.06, (a) the 3D phase portrait; (b) projected orbit on x-y plane; 
(c) projected orbit on x-z plane; (d) projected orbit on y-z plane. The Lorenz attractor corresponding to 
a = in the system ([3|, (e) the 3D phase portrait; (f) projected orbit on x-y plane; (g) projected orbit on 
x-z plane; (h) projected orbit on y-z plane. 



(a) 




(g) 12 (h) 



Figure 7: Attractor of system ([5| with r = 0.7, (a) the 3D phase portrait; (b) projected orbit on x-y plane; 
(c) projected orbit on x-z plane; (d) projected orbit on y-z plane. The Chen attractor corresponding to 
a = 1 in the system ([3|, (e) the 3D phase portrait; (f) projected orbit on x-y plane; (g) projected orbit on 
x-z plane; (h) projected orbit on y-z plane. 



5 Concluding remarks and future research 



5.1 Concluding remarks 

Regarding the few questions raised in the Introduction section, we now come out with some answers: 

1. In this paper, we have reported the finding of a structurahy simple yet dynamicahy complex one- 
parameter family of 3D autonomous systems with only two quadratic terms like the Lorenz system, 
with some key constant coefficients set as 1 or —1, leaving only one tunable real parameter in the linear 
part of the system. 

2. By tuning the singe parameter, this new family of systems can generate a chain of chaotic attr actors in 
a very similar gradual changing process like the one from Lorenz attractor to Chen attractor generated 
by the generalized Lorenz systems family and the one by the unified chaotic system. 

3. Unlike the typical Chen system, however, the new system that generates the Chen-like attractor does 
not satisfy the condition ai2a2i < in the linear part of the generalized Lorenz canonical form defined 
in [8]. 

4. Unlike the classical Lorenz system, the new system that can generates a similar Lorenz-like attractor 
does not satisfy the condition aiia22 > in the linear part of the generalzied Lorenz canonical form 
defined in [16]. 

5. This system is richer and more interesting than the unified chaotic system, in the sense that the new 
single parameter has a wider range to generate similar form of chaos. 

6. Similar to the generalized Lorenz systems family, the new system also has three equilibria, with one 
saddle and two saddle-foci, where the two saddle-foci have the same eigenvalues (one is negative real 
and two are complex conjugate with a positive real part). However, the in the new system the real 
part of the eigenvalues starts from zero, which can be used to define a starting point for the Lorenz-like 
systems family. In this sense, the new system has literally extended the unified chaotic system, even 
the generalized Lorenz systems family, to some extent. 

5.2 Future research 

Some related research issues may be further pursued in the near future: 

1. There may exist other Lorenz-type systems with the same nonlinear terms as the classical Lorenz 
system, which can generate a variety of chaotic attractors. 

2. There may even exist other types of nonlinear terms that satisfy the expected z-axis rotational sym- 
metry, for which the nonlinear term in the second equation must be like xz or yz^ and that in the third 
equation must be like xy^ or z^, or x^, or y'^ . So, if one wants to study 3D autonomous systems with 
two quadratic terms that can maintain the z-axis rotational symmetry, then there are totally 2x4 = 8 
possible combinations to consider. It is therefore interesting to ask if there would be other chaotic 
systems with two nonlinear terms different from those studied so far around the Lorenz system. 

3. In the new family of chaotic systems studied in this paper, the real part of the complex conjugate 
eigenvalues at system equilibria starts from zero for the case with r = 0.05, which can be defined as 
the starting point of the Lorenz-like systems family. This particularly interesting critical system with 
r = 0.05 alone deserves further investigations. 
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4. This paper has further revealed that some systems with different structures can hterahy generate 
similar-shaped chaotic attractors. Therefore, for 3D autonomous systems with two quadratic terms, 
the relation between the system algebraic structure and system chaotic dynamics is an important and 
interesting issue to be further revealed, understood and analyzed. 
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Table 1: Equilibria and eigenvalues of several typical systems. 



Systems 




Equations 




Equilibria 


Eigenvalues 


Unified system 
a = 


X = 10{y — x) 
y = 28x — y — xz 
z = -^z + xy 




(0,0,0) 
(±6 A ±6 A 27) 


-22.8277,-2.6667,11.8277 
-13.8546, 0.0940 zb 0.1945i 


Unified system 
a = 0.5 


y 


x=f{y-x) 
= ^x ^ ^y - xz 
z = -^z-\-xy 




(0,0,0) 
{±V68, ±V68, 24) 


-28.1696,-2.8333,19.1696 
-16.9593, 2.5630 zb 13.1857i 


Unified system 
a = 0.8 




X = 30{y — x) 

y= ^y-xz 

z = - j^z + xy 


(0,0,0) 
(±2v^,±2vf^,22.2) 


-30.000,-2.9333,22.200 
-17.9535, 3.6101 zb 14.3037i 


Unified system 
a = 1 


y 


X — 3b{y — x) 
= —7x + 2Sy — xz 
z = —3z + xy 




(0,0,0) 
(zb3x/7, zb3x/7,21) 


-30.8357, -3, 23.8359 
-18.4288, 4.2140 zb 14.8846i 


New system 
r 0.05 


y ■ 


X = —X — y 
= —X + O.Oby — xz 
z = —O.lz + xy 




(0,0,0) 
/42/20,T\/42/20,-1.05) 


-0.1,0.6544,-1.6044 
-1.05,0zb0.4472i 


New system 
r = 0.1 


y 


X = —X — y 
= —X + O.ly — xz 
z = —O.lz + xy 




(0,0,0) 

/n/io,TA/iT/io,-i.i) 


-0.1,0.6913,-1.5913 
-1.0162, 0.0081 zb0.4652i 


New system 
r = 0.5 


y 


X = —X — y 
= —X + 0.5y — xz 
z = — O.I2: + Xi/ 




(0,0,0) 
/I5/10,T^/i5/10,-1.5) 


-0.1,1,-1.5 
-0.8101, 0.1051 zb0.5994i 


New system 
r = 0.7 


y 


X = —X — 2/ 
= —X + 0.7?/ — xz 
z = — O.I2; + xy 




(0,0,0) 
/T7/10,T^/10,-1.7) 


-0.1,1.1624,-1.4624 
-0.7446, 0.1723 zb0.6534i 


New system 
r = 0.75 


y ■ 


X = —X — y 
= —X + 0.75?/ — xz 
z = — O.I2: + X?/ 




(0,0,0) 
/70/20,=F\/70/20,-1.75) 


-0.1,1.2038,-1.4538 
-0.7312, 0.1906 zb0.6651i 



Table 2: Numerical results for some values of the parameter r with initial values (1, 1, 1). 



Parameters 


Eigenvalues 


Lyapunov Exponents 


Fractal Dimensions 


r = 0.05 


Ai = -1.05 
A2,3 = ±0.4472i 


Li = 0.03467 

^2 =0 
L3 = -1.0843 


Dl = 2.0317 


r = 0.1 


Ai = -1.0162 
A2,3 = 0.0081 zb 0.4652i 


Li = 0.03058 

1.2 = 
L3 = -1.0301 


Dl = 2.029 


r = 0.5 


Ai = -0.8101 
A2,3 = 0.1051 =b0.5994i 


Li = 0.10531 

^2 = 

L3 = -0.70561 


Dl = 2.1497 


r = 0.7 


Ai = -0.7446 
A2,3 = -0.1723 =b0.6534i 


Li = 0.08953 

^2 = 

L3 = -0.48973 


Dl = 2.1832 


r = 0.75 


Ai = -0.7312 
A2,3 = 0.1906 =b0.6651i 


Li = 
L2 = -0.04626 
L3 = -3.0403 


Dl = 1.0063 
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